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A Study of Charmonium Systems across the Deconfinement Transition * 
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We present results from lattice studies of charmonium systems near the deconfinement transition temperature. 
On quenched isotropic lattices with lattice spacings between 0.02 and 0.05 fm, qq systems with quark masses 
close to the charm mass and with different spin-parity quantum numbers are studied in the temperature range 
0.9T C — 3T C . Results for temporal correlators of local operators, and the spectral functions constructed from them, 
are discussed. For the pseudoscalar and vector channels, the correlators are observed to change very little across 
the deconfinement transition, unlike in the case of the light quarks. 
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The behavior of charmonia across the decon- 
finement transition has been a subject of consid- 
erable interest ever since the breakthrough paper 
of Matsui and Satz jjj. Unlike light mesons, 
charmonia may exist as bound states even after 
the deconfinement transition. However, based on 
nonrclativistic arguments, Matsui and Satz pre- 
dicted that already at temperatures close to tran- 
sition, binding between quarks is reduced enough 
to dissolve J ftp, and proposed its suppression as 
a signal of deconfinement. Several later studies, 
based on potential model calculations, predicted 
a pattern of dissolution, with the higher excita- 
tions dissolving earlier, and J/ip dissolving at a 
temperature « 1.1 T c [|J. 

A more direct and reliable way to study the in- 
mcdium modifications of charmonia is to study 
on lattice the finite temperature imaginary time 
correlators 



G H (r,p,T) = (J ff (r,p)4(0, ~p)h 



(1) 



where Jh is the suitable mesonic operator, pro- 
jected on the state with spatial momentum p. For 
zero temperature studies, one often uses smeared 
operators to get a good overlap with the ground 
state. However, since we are interested in study- 
ing whether the bound states exist or not, using 



*This work was partially supported by DFG grant FOR 
339/1-2, and by GSI under the project BI-KAR. 
^Presented by S. Datta. 



smeared operators is not a good idea since it can 
mimick bound states j|. We therefore use point- 
to-point correlators, and study the temporal cor- 
relators, Eq. (|]), with Jh = cc, C75C, cj^c 
and C7 M 7 5 c to explore the 3 P (Xco), (Vc), 
3 Si (J/ip) and 3 Pi (xci) channels, respectively. 
For this study we use only quenched lattices, and 
used the nonperturbatively improved Q clover 
action for the quarks. Our lattice parameters are 
given in Table [I]. The lattice operators are con- 
nected to the continuum operators as J^ ont = 

z HflJ H ata ~ 3 , where /« = a/2k(1 + am ) is the 
quark renormalization factor for massive quark 
with tadpole improved bare quark mass afriQ = 
(1/2/c — 1/2k c )/u . Zh is obtained from tad- 
pole improved perturbation theory in the mass- 
less limit [|J. 



Table 1 

Lattice parameters. 



The lattice spacing is ob- 



taincd from the string tension. 


P 


a-^GeV) 


Size 


T/T c 


#conf 


6.499 


4.042 


48 3 x 16 


0.93 


50 


6.499 


4.042 


48 3 x 12 


1.25 


50 


6.499 


4.042 


48 3 x 10 


1.5 


46 


6.640 


4.860 


48 3 x 12 


1.5 


60 


7.192 


9.720 


48 3 x 12 


3.0 


90 



For our coarsest lattices, we have three temper- 
atures at the same lattice spacing. We use two 



k values, 0.1300 and 0.1234, which bracket the 
charm quark mass. We estimate the zero temper- 
ature masses of the corresponding mesons from 
the spatial correlators at 0.93T C . This is permis- 
sible since the spatial correlators correspond to 
zero-temperature correlators for a lattice with an 
extent of w 0.9 fm in one of the spatial directions, 
which is much larger than the typical size of a 
charmonium. (In fact, even for the much larger 
light mesons, the spatial correlators at this tem- 
perature have been found to reproduce the zero 
temperature masses quite well 0.) The masses 
we get are given in Table |[ By linear interpo- 
lation, we estimate that K c harm ~ 0.1280 for this 
lattice. 

Table 2 

Masses (in GeV) of the different representations 



Fig. |T] shows the temporal correlators for the 
different channels, for k = 0.1300. To highlight 
the difference from free quark propagation, we 
have divided the correlators by the corresponding 
free correlators, with same bare quark mass (in 
lattice units) as the tadpole improved bare quark 
mass of the interacting case. Fig. [l] shows that 
the scalar and axial vector channels are consider- 
ably modified on crossing T c . In the pseudoscalar 
and vector channels, on the other hand, there is 
very little change between 0.93 T c and 1.5 T c . 

Further insight can be obtained from looking at 
the spectral function at each temperature. The 
spectral function a(u>) is connected to the imagi- 
nary time correlator G(r) by the integral equation 



G(r) 



du>a{uj 



cosh(w(r - 1/2T)) 
sinh(w/2T) ' 



(2) 



(Here we consider only states projected to zero 
spatial momentum.) The severely undercon- 
strained problem of reconstruction of the spectral 
function from the information of the correlator 
at finite number of points can be handled using 
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Figure 1. Temporal correlators for the different 
channels (top: pseudoscalar and vector; bottom: 
scalar and axial vector) at j3 = 6.499, for k = 
0.1300. 

the maximum entropy method. At zero tempera- 
ture this method has been successfully applied Q| . 
At finite temperature, this method was also ap- 
plied in the light quark sector ||] . We follow here 
the methods of Ref. (^]. The spectral functions 
for the pseudoscalar and axial vector channels are 
shown in Fig. |2l While the axial vector peak be- 
low T c in Fig. ghas completely disappeared above 
T c , the pseudoscalar bound state peak persists at 
1.25 T c , with an essentially unchanged position, 
though somewhat broadened. The vector channel 
shows very similar behavior to the pseudoscalar 
one, with the peak persisting at 1.25 T c , while the 
scalar channel behavior is very similar to the ax- 
ial vector one. The pseudoscalar and vector peaks 
also seem to survive at 1.5 T c , though finer lat- 
tices are needed at these and higher temperatures 
to have a reasonably large number of points in 
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the temporal direction. The behavior of the dif- 
ferent channels is very similar in the k = 0.1234 
case, the scalar and axial vector channel peaks 
disappearing already at 1.25 T c while the pseu- 
doscalar and vector peaks surviving beyond this 
temperature. Since these two values of k bracket 
the charm, we conclude that for physical charm, 
while the IP states Xco and Xd are dissolved, the 
IS states J/ip and r\ c survive as bound states at 
temperatures up to 1.5 T c . 
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Figure 2. Spectral functions for the pseudoscalar 
(top) and axial vector (bottom) channels at 0.93 
and 1.25 T c . A detailed error analysis will be 
presented elsewhere. 

At the finer lattices listed in Table [l], we at 
present have only one temperature at each lat- 
tice spacing. So we need to estimate k for each 
case. We choose the k to give the same (tadpole- 
improved) bare quark mass, in physical units, as 
the K c harm obtained at j3 = 6.499. This way, we 
estimate k = 0.1290 at j3 = 6.64 (for our 1.5 T c 



run) and k = 0.13114 at /3 = 7.192 (the 3 T c run). 
The data at 1.5 T c supports the results of the 
coarser lattice at the same temperature. How- 
ever, since we can not estimate the systematic 
error involved in our estimate of n, more precise 
conclusion cannot be attained at present. We are 
now generating results at 1.1 T c for j3 — 6.64 and 
at 1.5 T c at (3 = 7.192. This will enable us to do a 
more detailed analysis of the results at these finer 
lattices. 

To summarize, we conduct a study of the be- 
havior of cc bound states across the deconfine- 
ment transition by looking at the suitable cor- 
relators and the spectral functions reconstructed 
from them. Our studies support the sequen- 
tial pattern for charmonium dissolution obtained 
from potential model studies, where the broader 
bound states (the scalar and axial vector chan- 
nels) dissolve before the pseudoscalar and vector 
channels. The pseudoscalar and vector channels 
are seen to survive as bound states still at 1.25 
T c , and probably dissolve after 1.5 T c . Results of 
direct lattice studies of charmonia systems have 
also been presented by M. Asakawa and K. No- 
mura in this conference ||]. Our conclusions are 
consistent with theirs. 
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